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ABSTRACT 


It is now understood that deep convection in the ocean plays a dominant role in 
determining the climate of the world's oceans. Recent theoretical advances in explaining 
oceanic convection need to be tested by real observations. Satellite observations of deep 
convection regions may be a promising new tool in studying this phenomenon. This thesis 
examines deep convection events in two ways: 

© To assess the characteristic elements of a deep convection event using two 

different prediction models based upon the turbulent kinetic energy budget. 

© To attempt to observe deep convection phenomena signals in altimeter data. 

In 1987, a deep convection event was observed in the northwestern Mediterranean 
sea (Schott and Leaman 1991). These data, combined with GEOSAT altimeter data, were 
used to verify the Kraus and Tumer and the Naval Postgraduate School mixed layer model 
predictions of the time evolution of temperature, salinity and mixed layer depth. Both 
models predicted final values similar to the observations, but model tuning was required to 
reproduce the observed rapid mixed-layer deepening. The interpolated altimeter field does 
not allow identification with confidence of the Mediterranean convection area. However, a 
locally persistent feature and the mean winter sea surface topography field agree with in- 
situ observations and do provide some indication about where and when the convection 


process occurs. 
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I. INTRODUCTION 


It is now recognized that deep convection to depths exceeding about one kilometer 
occurs in very few areas in the world's oceans. Locations of deep convection are mostly 


limited to polar regions such as the Labrador, Weddell, and Greenland seas. However, 


deep convection also occurs in the Northwestern Mediterranean sea. The process of deep 


convection is very important in deep-water formation and thus plays a fundamental role in 
determining the climate of the world's oceans. 

Deep convection in the open sea has been studied for more than twenty years. Ob- 
servations and models are now available. These models have helped explain the phenome- 
non on a global scale, but not the small scale. Recent theoretical advances have been made 
(Garwood 1991; Garwood et al. 1994; Paluszkiewicz et al. 1994) and need to be tested by 
real observations. Satellite observations of deep convection regions have rarely been ana- 
lyzed and the large variety of remote sensing instruments may open promising new ways 
to study this phenomenon. 

The 1987 Mediterranean convection experiment (Schott and Leaman 1991; Lea- 
man and Schott 1991) documented a deep convection event in the Golfe du Lion. Data 
were recorded between January 25 and March 6. Altimeter data from GEOSAT for the 
same period and location, but independent from the experiment, were also used in this 
study. 

The purpose of this thesis is examine a Mediterranean deep convection event in two 
ways: 

© To assess the characteristic elements of a deep convection event using two one- 

dimensional models based upon the turbulent kinetic energy (TKE) budget. 

© To look at deep convection phenomena using remote sensing instruments. 

The Kraus and Turner (1967) model and the Naval Postgraduate School model 
(Garwood 1977; 1991) with realistic forcing and initial conditions, will be used to predict 
the evolution of temperature, salinity, TKE and mixed layer depth with time and to 





compare them with in-situ profiles. An analysis of the satellite data will be conducted to 
identify sea level changes in time and space, thus potentially being able to remotely identify 


the deep convection process. 


I. THE CONVECTION PROCESS IN THE MEDITERRANEAN SEA 
A. A THREE-PHASE PROCESS 

A few areas in the world ocean have been identified as locations of winter deep 
convection (Killworth 1976). Along with the Labrador, Weddell and Greenland seas, the 
northwestern Mediterranean Sea, especially the Golfe du Lion south of France, is also a 
convection-favorable region. Due to easy access and often predictable triggering factors, 
the deep convection and the formation of deep water in the Golfe du Lion have been 
the subject of a few observational studies and modeling. From 1969 (Medoc group 1970) 
to the 1987 experiment (Schott and Leaman 1991, Leaman and Schott 1991), deep con- 
vection has been observed many times, each time improving the understanding of the "'pre- 


conditioning" phase, the "violent mixing" phase and the "sinking and spreading" phase 


(Medoc group 1970) of deep convection (Fig. 2). 























Figure 1. Bathymetry of the Golfe du Lion, in Northwestern Mediterranean Sea and the 
Locations of the CTD and PEGASUS Stations and the ADCP Moorings during the 1987 
Mediterranean Convection Experiment From Leaman and Schott 1991. 





The pre-conditioning phase is characterized by three favorable factors (Fig.1; 
Schott and Leaman 1991). Firstly, the wintertime cyclonic circulation induces a doming of 
the isopycnals in its center causing the typical Mediterranean three-layer stratification to 
have a noticeably weaker stratification below the surface layer (Swallow and Gaston 
1973). The second factor is the conelike topography of the Rhone fan (42°N-5°E), which 
may induce cyclonic circulation with a predictable position for the center of the dome each 
winter (Hogg 1973). The last factor is the incidence of two cold and dry winds, originat- 
ing from the continent and directed toward the center of the dome: Mistral from the 
Rhone Valley and Tramontane from the Pyrenees. Cooling and evaporation are increased 
by these winds which act to destabilize the oceanic upper layer (phase 1 in Fig. 2). Once 
this convective instability occurs, the preconditioning phase is followed by vertical convec- 
tion associated with large downward vertical velocities (Fig. 3) which penetrate to deep 
levels since the pre-existing stratification below the surface is weak. During the last phase, 


the new deep water spreads away from the formation area. 
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Figure 2. Schematic of a Deep Convection Process: Evolution of the Density Profile. 


B. DEEP CONVECTION MECHANISMS AND MODELS 

The mechanisms responsible for the establishment of the mixed water column have 
been the subject of numerous investigations. The mechanisms proposed include gravity- 
inertal waves (Gascard 1973), convective process with baroclinic adjustment (Gascard 
1978), comparison to atmospheric thermals (Vaughan 1993) or parcel instability mecha- 
nism (Garwood 1991), taking into account "thermobaricity" (increase of the thermal ex- 
pansion coefficient with pressure). A few models have been developed to predict what the 
MEDOC group in 1970 characterized as a homogeneous column of water created by ver- 
tical mixing due to strong surface cooling. Killworth (1976) used a two-dimensional 
model where the violent mixing phase was explained by one dimensional non-penetrative 
convection coupled with geostrophic adjustment. A multilevel numerical model (Madec et 
al. 1991) of the Laboratoire d'Oceanographie Dynamique et de Climatologie (LODYC) 
used the mechanism proposed by Gascard (1978). 


Cc THE RESULTS OF THE MEDOC EXPERIMENT 

Medoc 69 (Medoc group 1970) and the experiment of the winter of 1987 (Schott 
and Leaman 1991, Leaman and Schott 1991) provide the most extensive observations of 
convection in the Northwestern Mediterranean Sea. In 1969, a large (20 km x 50 km) and 
deep (1500-2000 m) well-mixed region was observed over the Rhone fan. The vertical ve- 
locities, measured by neutrally buoyant rotating floats, showed upward and downward ve- 
locities O(2-3 cm/s). The existence of rotating plumes O(6 km diameter) was also 
deduced (Voorhis and Webb 1970). 

Using Acoustic Doppler Current Profilers (ADCP), Vector Averaging Current 
Meters (VACM), Conductivity Temperature and Depth instruments (CTD) on site and 
coastal meteorological measurements, the 1987 experiment collected numerous data. 
ADCP measures currents in three dimensions over a fairly large range with an accuracy 


and time resolution allowing measurements of vertical velocities during convective events. 


An array of three moorings with upward-looking ADCP on top was deployed in the Golfe 





du Lion between January 24 and March 6. Unfortunately, a strong Mistral had already oc- 
curred in early January and triggered the convective process. A second convective event 


took place between February 12-22 and was fully recorded. 
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Figure 3. Time Series of Vertical Component of the Current at 375 meters in February 
1987. 


The time series of the vertical current in Fig. 3 shows large downward velocities 
up to 10 cm/s associated with the second Mistral period in mid February 1987. There was 
almost no change in the velocity records with depth, confirming that the violent mixing 
phase is very swift and occurs throughout all the water column. The variability of the ver- 
tical kinetic energy observed during the calm and the windy periods tends to refute the in- 
ternal wave, mechanical wind or wave forcing theories. This variability is associated with 
convection following surface cooling (Schott and Leaman 1991). Based on the assump- 
tions that (1) the downward velocity events are quasi-stationary at the time scale of the 


observations (one hour) and (2) that the background flow is homogeneous, a diameter of 


1 kilometer can be deduced for the convection cells (Schott and Leaman 1991). They are 
embedded in the convective area itself O(100 km) from which eddies O(5 km) detached 
(Gascard 1978). Other interesting results concern the surface circulation and the move- 
ment of the lighter surface water between Mistral events. A correlation was also found be- 
tween deep water temperature changes and the magnitude of surface cooling. Finally, 


gravitationally unstable surface layers were observed during periods of heightened atmos- 


pheric forcing. 





Il. ASSESSMENT OF CHARACTERISTIC DEEP CONVECTION ELEMENTS 
USING TWO SIMPLE MODELS 


In this chapter the characteristic elements of a deep convection process are as- 
sessed. The Kraus and Turner model and the Naval Postgraduate School model will be 
used to evaluate the deepening and the evolution of temperature, salinity and density of 
the mixed layer during that phenomenon. Both models will be initiated with real data and 
the results of the runs will be compared to the observed profiles. The Kraus and Turner 
(KT) model and the Naval Postgraduate School (NPS) model are both one-dimensional 
models based on the vertically integrated turbulent kinetic energy (TKE) equation. They 
differ by the way some of the terms of the TKE budget are modeled or neglected, the KT 


model being the simpler of the two. 


A. THE GOVERNING SET OF EQUATIONS 


1. Zones of the water column 





‘emperature 
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Figure 4, Typical Temperature Profile in a Water Column. 


The following equations refer to a typical water column which can be divided into 
three regions as seen on Fig. 4: 

© From the surface to depth h, the well mixed layer is where the temperature and 
salinity are constant due to strong turbulent mixing. 

© The interface region between h and h+8, is a region of sharp temperature and sa- 
linity gradient also called entrainment region. 

© The turbulence free or stratified zone is where the temperature decreases with 
depth. 

2. The temperature and salinity equations 

The equation of conservation of mean temperature with the assumption of hori- 


zontal homogeneity (vertical gradients much larger than the horizontal gradients), yields: 


a wear, a a 
yi + WS + 5 Tw) =0 (3.1) 

We first integrate Eq. 3.1 with respect to z over the mixed layer depth (MLD) h. 
There is no advective term because T is constant in the mixed layer, and we use the as- 
sumption that the entrainment zone thickness 5 << h. Knowing that the surface heat flux 


can be rewritten using the total heat flux Qo, the density p and a heat coefficient Cp: 


WP(0)= 22 
= ()= cep (3.2) 


the combination of Eq. 3.1 vertically integrated and Eq. 3.2 gives a prediction equation for 


the mean temperature in the mixed layer: 


F_ Qo ws 


of pCph sh 6.3) 


In a similar manner, we have an equation for the mean mixed layer salinity: 


oS _ S(e-p)_wS'4 

ok h Ge) 
where(e-p) is the net evaporation minus precipitation. 

The rate of mixed layer deepening with time is obtained by integrating Eq. 3.1 


across the interface. Since the interface is moving, we have to use the Liebniz rule, 


aa) 
or 


db(t) 


fe Lae = 2 I" Face Mz= an GO -re= 625) 


a) Ot Ot Jal 


Using F = T and limits such that a(t) = h(t) and b(t) = h(t)+8, this integration gives 
the "jump condition" which relates the entrainment velocity to the temperature jump at the 


base of the mixed layer and the heat flux across the interface, 


(3.6) 





This means that h can be changed by entrainment and upwelling or downwelling. 
A jump condition can be also written for the salinity : 


_ ws 
w= eS (3.7) 


We shall assume here that there is no up/downwelling so that the change with time 
of the mixed layer depth is only due to entrainment. 

The combination of Eq. 3.3, Eq. 3.4, Eq 3.6 and Eq. 3.7 gives new prediction 
equations for the mean mixed layer temperature (3.8) and salinity (3.9). 


aT__Qo_ ATah 
ot p.Cp.h = h Or G8) 





OS _ S(e-p) _ AS.we (3.9) 


This equation for the mean sea surface temperature simply indicates that both sur- 
face heating/cooling and entrainment mixing will change the mean temperature in the 
mixed layer. 

3. The turbulent kinetic equation 

Our system is now composed of two equations (Eq. 3.8 and Eq. 3.9) and three un- 
knowns: 7,5 and h or we = %. To solve it, we need another equation which will be pro- 
vided by the turbulent kinetic energy (TKE) equation: 


LENE i aw Ws 
SER) = HW E-VWE) + (Boe) ~ Sw A + Fe B.10) 
storage shear buoyancy turbulent transport dissipation 


production flux 
There is assumed to be no turbulence below the base of the mixed layer. 


B. |THE KRAUS AND TURNER CLOSURE (1967) 

Kraus and Turmer used the TKE budget with the following assumptions to derive 
the third equation: 

© The pressure redistribution in the turbulent transport term is neglected. 

© The dissipation is included in the shear production and buoyancy flux through 

the tuning coefficients (c, and c,). 
© A steady state is assumed. 
The equation of state is used to express the buoyancy flux in terms of heat and sa- 


linity fluxes. 


P= Poll —a(T-T)+B(S—So)] (3.11) 


12 


With these assumptions, the vertical integration of the TKE equation reduces to : 


ct cafe 22+ Beste) (aga T-BeAs3t=0 (3.12) 


h 

‘The integrated transport of TKE averaged over the mixed layer is zero because 
there is no external source or sink for turbulence. This equation gives the entrainment ve- 
locity or rate of deepening in terms of wind stress, heat and salinity fluxes and temperature 
and salinity jumps at the base of the mixed layer. 

For our purpose we will assume that the contribution of the salinity term is much 
smaller than the heat flux term. This can be proven by a quick computation. The excess of 
evaporation over precipitation in Mediterranean sea can be estimated to | meter per year 
(Pickard et al. 1990). Thus the salinity term in Eq. 3.12 has an order of magnitude of 
400 m’/h> whereas the heat flux term has a mean value of 4000 m’/h’ during the convection 
period. This ratio of 10 is enough to neglect the salinity contribution. 

Equations 3.8, 3.9 and 3.12 are now used in a Matlab program to evaluate the 


evolution of the depth, temperature and salinity of the mixed layer with time. 


C. | THE KRAUS AND TURNER SIMULATION 

1. Initialization of the simulation 

The Kraus and Tumer model is simple and easy to program. The input variables 
are the initial potential temperature and salinity profiles, the total heat flux, the friction ve- 
locity and the tuning coefficients. 

Two series of simulations will be examined. First, both coefficients will be set to 
unity to give an equal contribution by the wind and heat flux. We will examine cases with 
both constant and variable atmospheric forcing. Then under a variable forcing, we will ad- 
just the coefficients to try to reproduce the Mediterranean convection event. 

The program runs for ten days with digitalized observations for initial temperature 


and salinity profiles (Fig. 5) and real forcing (Fig. 6). Notice that the strong unstable 
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temperature profile is balanced by a strong stable salinity profile to produce a typical 
Mediterranean winter profile (stable-neutral at the surface / strongly stable for the interme- 


diate layer / neutral for the deep layer). 
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Figure 5. Digitalized Temperature, Salinity and Density Profiles at the Beginning of the 


Deep Convection Event 


This ten day period corresponds to the Mediterranean deep convection event ob- 
served in February 1987 (Schott and Leaman 1991; Leaman and Schott 1991). The period 
was delimitated by examining the period of strong atmospheric forcing (February 12 - 22) 
that corresponds to the period of maximum downward velocity as well as homogenization 


of the water column. 
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Figure 6. Total Heat Flux and Windstress vs. Time during the Deep Convection Event. 
After Leaman and Schott 1991. 


The Matlab function ode45 was used to integrate the system of ordinary differen- 
tial equations using 4th and Sth order Runge-Kutta formula. 
The results of the different runs are discussed in the following section, the plots are 


presented at the end of the chapter. 
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2. Constant coefficients simulation 

a. Constant atmospheric forcing 

The constant values used as input for the atmospheric forcing are the 
means of the total heat flux and friction velocity over the ten-day convective period. The 
mixed layer deepens to 950 meters after ten days as shown on Fig. 11. There is an inflec- 
tion point around hour 75 that indicates an extrema in the slope of the curve, thus a 
change in the rate of deepening. Since the forcing is constant, the change is linked to the 
change in temperature and salinity profiles at 250 meters. 

The mixed layer temperature shows an early peak that can be explained by 
the mixing of the intermediate warm water with the colder surface water. This mixing also 
takes 75 hours to complete. When the mixed layer is then homogenized, the situation be- 
comes similar to the classical two-layer mixing case in which warm water overlays cold 
water. The temperature decreases uniformly to 12.79°C. Meanwhile, the mixed layer salin- 
ity steadily increases, which is normal when mixing the initial profile. The density profile 
shows no significant change during the period. 

b. Variable atmospheric forcing 

Allowing the wind stress and total heat flux to change with time (Fig. 12) 
produces its most visible effect on the mixed layer depth. The depth reached after ten days 
is now 1150 meters. The curves are not as smoothed as in the previous plots and reflect 
somewhat the variations in the forcing. However, the global shape of all curves is 
preserved. 

During the first 50 hours, the forcing is weak, thus the curves show a slow 
departure from the initial conditions. Since the deepening is slower, the mixing with the 
warmer and saltier intermediate water takes more time with peak temperature at 80 hours 
instead of 20 hours. On the other hand, the values of more intense surface forcing found 
around hour 200 induce a rapid increase in the mixed layer depth. At that time, the tem- 
perature and salinity have already reached their final values and show no further 


evolution. 
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c. Comparison to in-situ profiles 

The in-situ deep convection temperature and salinity profiles of Fig. 7 
(Leaman and Schott 1991) show a rapid homogenization of the layer. Within a few hours 
to half a day, the stratification has totally disappeared, and the profiles are simply straight 
lines. This also underlines the importance of the preconditionning phase. The water col- 
umn is so well prepared that the weak forcing at the beginning of the period (February 13) 


can trigger the deep convection process. 
In our simulation, the KT model deepens the layer to the bottom in more 


that ten days. The maximum depth reached was 1200 meters with a bottom depth of 2500 
meters. However, the final mixed layer temperature and salinity (12.79°C and 38.445) 
agree fairly well with the temperature and salinity found in the water column immediately 
after the convection event (12.81°C and 38.45; Fig. 7). 
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Figure 7. In-situ Potential Temperature, Salinity and Sigma Theta Profiles before and after 


a Deep Convection Event. From Leaman and Schott 1991. 





The temperature recorded at 537 meters during the convection has been 
averaged over six hours and is presented in Fig. 8, along with the original recording. This 
measurement can be assimilated to an Eulerian observation. The time series goes further in 
time than the profiles of Fig. 7. Some gaps in the data recording are present but small time 
scale variations are evident while the temperature decreases slightly. Although the varia- 
tions are very small and must be analyzed carefully, this indicate that the convection proc- 
ess is not a simple event in time but a sum of short time scale events and/or that even 


when the water column seems homogenized, there is still turbulent mixing occurring. 
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Figure 8. Potential Temperature at 537 Meters during the Convection Event 


On the contrary, the simulated temperature is the mean temperature above 
the interface, the time series represents the temperature evolution as the interface deepens. 
It is a Lagrangian measurement. The KT physics do not account for subsequent cooling 


that would modify the upper layer and create a new colder and sinking water mass. 
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The Fig. 9 is what a sensor at 530 meters would see, using the information 


we know about temperature and depth from the KT simulation. 





‘Simulated potential KT temperature at 530 meters 
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Figure 9, Simulated Temperature at 530 meters during the Convection Event 


Unfortunately, there is no possible comparison. The simulated temperature 
evolves only after the interface has reached the level of 530 meters, almost at the end of 
the convection period. This revels some of the limitations of this model; successive cool- 
ings and the rapidity of the process cannot be modeled. 

This model extends the process of deep convection over a longer period 
than found in reality. The final results are consistent with the real case and are reasonable 
considering the simplicity of the model. 

3. Variable coefficients simulation 
Simples cases such as cl=0 (heat flux forcing only) or c2=0 (wind stress forcing 
only) will allow the assessment of the respective roles of the wind stress and heat flux in 


the deepening of the mixed layer. Then we will double the coefficients one at a time, and 


finally all together. Since were are using a relatively simple model, our aim is to have the 
mixed layer depth reach the bottom by the end of the period of strong surface forcing. 

a. Heat flux forcing alone (c1=0) 

The deepening due to surface cooling alone leads to a mixed layer depth of 
650 meters (dot dashed plots in Fig. 13). The first two days of weak heat flux are evident 
in all the plots: the rate of deepening is zero and there is no evolution in either the tem- 
perature or the salinity. When the heat flux increases it triggers the mixing, producing pro- 
files similar to those found in the previous cases. The changes in the deepening rate 
correspond to the variations in the heat flux forcing at hours 100 and 175. Beside the 
mixed layer depth, all the final values are consistent with the observed ones. 

b. Wind forcing alone (c2=0) 

In the KT simulation, the wind (solid plots in Fig. 13) doesn't seem to be 
sufficiently powerful to mix to a deep level. The mixed layer depth reaches 275 meters. 
The wind was significant during that period (10 to 20 m/s), and its variations can easily be 
seen in all the curves. Since the rate of the mixing is slowed, the peak and final tempera- 
tures are both slightly higher than for the previous runs (12.915°C and 12.888°C, respec- 
tively). For the same reason, the salinity and density profiles show lower final values. 

Therefore, in the variable atmospheric forcing, the heat flux forcing is the 
dominant mechanism responsible for the mixing. 

c. Adjustment of the tuning coefficients 

The adjustment is based on the variation of the mixed layer depth in time. 
The temperature, salinity and density profiles are similar for the three following runs; they 
show exactly the same shape as for the constant coefficients section and reach the ob- 
served values after ten days. 

© Experiment #1: cl=1 / c2=2: We have already shown that the heat flux was the 
dominant process of the mixing. Thus by increasing its coefficient, the bottom is 


reached sooner ( dotted plots in Fig. 14). 


20 


© Experiment #2: cl=2 / c2=1: The result is not very different from the original 

variable atmospheric forcing case with both coefficients set to unity. Again, this 

tells us that the influence of the wind stress is weak (dot dashed plots in Fig. 14), 

© Experiment #3: cl=2 / c2=2: It is the combination of the two previous experi- 
ments. The bottom is reached at time 215 hours (solid plots in Fig 14). 

This set of tuning experiments proves that the bottom can be reached with 

a KT simulation, This is easily done by doubling the c2 coefficient, increasing the heat flux 

influence. This may also signify that in a deep convection process the heat flux largely 


dominates the wind stress, however strong it is. 


D. THE NAVAL POSTGRADUATE SCHOOL CLOSURE 
The Naval Postgraduate School mixed layer model (Garwood 1977) is a one- 
dimensional model using a bulk second-order turbulence closure system for closure using 
TKE equation. This model differs from the Kraus and Turner model in five ways: 
© Three equations now describe the TKE budget (one for each component u” , v? 
and w”). 
© The dissipation is equally distributed into the three parts of the TKE equation 
(D/3 for each equation; this is a good approximation for turbulent flow). 
© The pressure redistribution is modeled (Ry, Ry and R,). 
© An additional shear term appears when we integrate the shear production term 
of the TKE across the interface (neglected in KT). 
* The horizontal momentum equations are now needed to compute the additional 
shear production term. 
1. The NPS mixed layer model equations 
The three vertically integrated turbulent kinetic equations are : 


hu 
or ri 


7 
= 2M SPX + auywe+R.—2 





EPs aviwetR,—2 G.15) 





cp * Bahsle— ‘p)- gh(aAT— BAS +R. ~2 (3.16) 


The right hand sides of Eq. 3.14 and 3.15 are the sum of several terms. The first 
term is the windstress part of the shear production where (t,, t,) are the components of 
the wind stress, p the density and u* the friction velocity. It is followed by an additional 
shear production term associated with the jump in velocity (AU,AV) between the well- 
mixed layer and the entrainment zone. R,,,., and D are pressure redistribution and dissipa- 
tion terms respectively. In Eq. 3.16, the entrainment buoyancy flux depends on the coeffi- 
cients @ and for thermal and salinity expansion, the temperature and salinity jumps AT 
and AS between the well-mixed layer and the entrainment zone. 

‘The mixed layer temperature and salinity are given by Eq. 3.9 and Eq. 3.10. 

The horizontal momentum equations are : 


we A = flavyt (3.17) 


oy) =f) +2 (3.18) 


a = 


The mixed layer depth equation is given by ors 





where w., is the mean vertical velocity at the base of the mixed layer. 

The system is now composed of 8 equations and 9 unknowns (T, S, h, Ry, Ry, R, 
D, U and V). To compute the model, we will express some of the unknowns in term of to- 
tal turbulent kinetic energy (E) and use the TKE budget at the interface to get the ninth 


equation. 
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2. Special functions and constants 

The friction velocity is defined by u* = (EE)? +(By"}"* (3.20) 
‘The net downward surface heat flux Q.=Qyi.- Quer natinion" Qinest” Quensistes 
The total TKE is : E=u? + V7 + Ww? 


In cases of usual mixed layer deepening, the dissipation € is parametrized by 
£ where the dissipation scale tp is a fraction (M1) of the time scale for one eddy rota- 


tion: ¥* When mixed layers are very deep, the Rossby number is O(1) meaning that the 
E 


planetary rotation has some influence on the water column motion. This can be taken into 





account by adding a fraction (M5) of the rotational time scale tg = re This yields a verti- 
cally integrated dissipation term: 
0 = = 
D=J? ede=2M:E +MSphE (3.21) 


This says that dissipation will adjust to the supply of total kinetic energy E and is 
influenced by rotation. 


The components of the pressure redistribution are modeled using a "return to en- 


tropy" theory: 
Re= ame -—W)E Ry= ama(E—WE R= ame -WE (3.22) 


‘The pressure redistribution can be viewed as an energy exchange between the three 
components of the TKE. 
The key point in this model is the entrainment velocity w,. An equation for w, 


comes from the TKE budget at the interface. Since we are in a region where there is 
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negligible turbulence, the shear production and dissipation terms have sufficiently small 


contribution that they can be included in the remaining terms. 


‘The unsteadiness term x can be approximated by wel orw = 


The buoyancy flux term is now a big loss term and is rewritten using a fraction 


(M4) of the temperature and salinity jumps at the base of the mixed layer: 


= 
ree -g = M4lagT’W — BSW] = M4[(—-cgAT + BgAS)]we (3.23) 





Finally, the turbulent transport term is approximated by z with tw = 





y 
termined by the quantity of energy from Ee carried by w', Solving for w, , we have: 


wee wees (3.24) 

Some asymptotic conditions, such as no stratification or free convection, can be 
applied to the model for basic verification. The results show that this entrainment model 
agrees with the extrem cases. 

3. Role of each coefficients 

The M coefficients introduced in these equations will be used later to tune the 
mixed layer model for application to deep convection. M1 and MS quantify the disssipa- 
tion, respectively for the TKE and rotation contribution. M2 is the pressure redistribution 
coefficient, which only appears in the individual vertically integrated TKE equations. It 
quantifies the energy exchange within the TKE budget. M3 quantifies the role of the shear 
production and therefore the windstress. M4 weights the buoyancy flux and indirectly the 


entrainment velocity term. 
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E. THE NPS SIMULATION 

In that series of simulations, we will first set the references by looking at both con- 
stant and variable atmospheric forcing under normal conditions (specified later). Then we 
will briefly estimate the influence of the Coriolis parameter, and finally assess the role of 
each 'M' coefficient to try to predict the observed states. 


The initialization is done with the same initial conditions as the KT simulation (ini- 


tial temperature and salinity profiles of Fig. 5) and the program runs for the same ten-day 


period. 
1, Simulation under normal conditions 
The normal conditions refer here to a Coriolis parameter taken at the mean latitude 

of the convection area and the M coefficients determined for a usual mixed layer deepen- 
ing (non convection type). These coefficients are set to M1=1, M2=0.5, M3=6, M4=1 and 
M5=1 (Garwood 1977). 

a. Constant atmospheric forcing 

The constant values used as input for the atmospheric forcing are the ten- 
day heat flux and friction velocity means. Besides deepening to 450 meters (Fig. 15), the 
smoothed curves are similar to the simulation using the KT model. The rate of deepening 
is fairly constant over the period; the peak temperature is reached after one day, and the 
final values for temperature, salinity and density are again close to the observed ones but 
slightly lower than with the KT simulation. 

b. Variable atmospheric forcing 

The shape of all curves (Fig. 16) reflects the variations in the atmospheric 
forcing, being flat when the forcing is small and steeper when the forcing increases. This 
results in delaying the temperature peak because of the delay in the surface forcing. The 
final mixed layer depth is 500 meters and all the other final values are the same as for the 
constant forcing simulation. 

As in the previous section, the comparison of the simulated temperature at 


530 meters and the in-situ recording (Fig. 8) will not be conclusive. Since the mixing 





reaches only 500 meters, the plot would simply show a flat curve at 12.8°C for the simu- 
lated temperature at 530 meters! 

2. Role of the Coriolis parameter 

Using a variable atmospheric forcing and the same values for the M coefficients, 
let's try to see the influence of rotation.The most visible effect of rotation is to decrease 
the final mixed layer depth as the Coriolis parameter increases (Fig. 17). for f=0, the final 
depth is 1400 meters or more than twice the value of f > 0. With and without rotation, the 
deepening rate is the same up to 250 meters or five days after the beginning of the convec- 
tion, but it dramatically increases after when there is no rotation. The temperature and sa- 
linity jumps are above 250 meters. Therefore once the mixed layer depth reaches 300 
meters, all the mixing is completed. For that reason, the final temperature, salinity and 
density values remain the same as for the previous cases. 

An explanation of this interesting result can be found in Garwood (1977) who de- 
scribes the general solution of the NPS model in terms of dimensionless parameters. He 
stated that the nondimensional mixed layer depth Z’ is directly proportional to the inverse 
boundary-layer Rossby number and thus proportional to the Coriolis parameter. The 
boundary-layer Rossby number (Ro) is the ratio of the friction velocity u’ to the Coriolis 
parameter. The rate of entrainment decreases with an increasing Z* or an increase in layer 
stability. Then it takes more time for the mixed layer to deepen at mid-latitudes than at the 
equator. For the same time period the mixed layer depth will be larger at the equator than 
at mid-latitudes. We can also state that in deep convection, the overturning creates large 
eddies which can be influenced by rotation. Thus, when there is no rotation, the deepening 
is more important. 

3. The influence of the M coefficients 

a. Search for the best coefficients 
All five coefficients influence the rate of mixed layer deepening. The goal is 
to determine the correct magnitude or weight of each coefficient in order to succesfully 


simulate a deep convection event. A simple doubling or halving of each coefficient one at 
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a time already shows significant changes in the mixed layer depth and temperature, Figure 
18 shows the variation of these two quantities with respect to the variations of each coeffi- 


cient. The solid curve serves as a reference; it is taken from the previous run which was 


done for a variable forcing, a Coriolis parameter taken at the Mediterranean latitude and 


the M coefficients at M1=1, M2=0.5, M3=6, M4=1 MS5=1. 

All coefficients seem to have a important effect on the deepening rate. The 
mixed layer depth increases significantly for a decrease in M1, M4 and MS and an increase 
in M2 and M3. All the temperature curves have the same shape and final value. 

Is there a reasonable explanation for the way the coefficients must be 
changed? The dissipation is linked to the TKE by M1 and to the Coriolis parameter by 
MS. Since the convection occurs to deep levels, the TKE in the water column is very large 
and the rotation may have a large influence on the column. Decreasing M1 and MS is 
needed to compensate for the large quantity of TKE produced and to reduce the influence 
of rotation. From the Kraus and Turner simulation, it was found that the buoyant produc- 
tion was a more efficient source of mixing than the shear production. Therefore, we can 
conclude that (1) since at great depth the wind is no longer felt, M3 must increase to keep 
a relative influence of the wind; (2) consequently M3 has to increase to balance the terms 
of the TKE equations; and (3) M4 has to decrease to reduce the influence of the efficient 
buoyancy flux. Decreasing M4 increases the entrainment velocity which is the case during 
the convection events. 

b. ANPS simulation with deep convection coefficients 

An investigation of the best coefficients is now possible by each time dou- 
bling M2 and M3 and halving M1, M4 and MS. Other choices for the coefficients are pos- 
sible, but this one respects the initial ratio between M1, M2, M3, M4 and MS and favors 
the rapidity of the process. 

Plotted on Fig 19, the increase in the coefficients shows a MLD reached 
sooner, thus an earlier mixing with the warmer and saltier intermediate water. We also no- 


tice that in the case of deep convection coefficients, the mixed layer depth responds to the 





significant forcing occurring after the second day whereas the temperature, salinity and 
density are mixed within the first few hours. This is consistent with the definition of vio- 
lent mixing in the second chapter. The coefficients kept for the NPS deep convection case 
are M1=0.125, M2=4, M3=48, M4=0.125 and M5=0.125. 

c. Validity of the coefficients found 

There are a few ways to assess the validity of the coefficients. First by 
looking at the time series of the TKE provided by the model for the 2500 meters water 
column and comparing it with the mean in-situ TKE at 450 meters. As the coefficients 
evolve, Fig. 20 shows an increase in the mean value of modeled TKE as well as a better 
fit for the in-situ TKE. The increase in TKE must be compensated by a decrease of M1. 

Then going back to the way Garwood (1977) found the values for its coef- 
ficients, we can try to apply the same method in our deep convection case. The easiest co- 
200) where u’ is the friction velocity and 8U(0) the 
magnitude of the mean velocity associated with the mean shear. The computation using in- 
situ data for U and V gives M3 O(10). 


The simulated temperature at 530 meters with deep convection coefficients 


efficient to compute is M3. M3 = 


(Fig. 10) shows a better agreement with the observed result. 





Simulated potential NPS temperature at 530 meters 








700 750 200 250 
‘Time In hours / Convection period i} 








Figure 10. Reconstitued Temperature at 530 Meters with Deep Convection Coefficients. 
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After the interface has reached the level 530 meters, the reconstitued tem- 


perature shows a decrease, similar in magnitude to the observed one, but much quicker in 


time. This only proves that there was some improvement compared to the Kraus and 
Tuner method, but we still need a better fit of the coefficients. 

Using in-situ data and a computation of the TKE budget from Eq. 3.10, we 
can estimate the mean dissipation in the case of deep convection, providing values for M1 
and MS using Eq. 3.21. This yields M1 and M5 O(0.1) which is consistent with the previ- 


ous estimation. 


F. SUMMARY OF THE SIMULATIONS 





Forcing | Coriolis | Kraus & Turner Naval Postgraduate 
type Parameter model School model 





f=0 MLD 950 meters MLD 1250 meters 


Constant 
forcing 





MLD 450 meters 
Usual 





coefficients MLD 1150 meters MLD 1400 meters 








MLD 500 meters 











Optimized | Variable MLD 2500 meters MLD 2500 meters 
|coetficients| forcing | after 215 hours after 75 hours | 

| 1 i | 
Table 1. Summary of the Mixed Layer Depth Reached using the Kraus and Tumer Model 














and the Naval Postgraduate Schoo! Model. 





G. PARTIAL CONCLUSION AND RECOMMENDATIONS 

The two simple and one-dimension models were used to evaluate the long term 
evolution of the mixed layer characteristics. The final computed values are almost the 
same as the observed values. However, with the tuning coefficients set for a normal mixed 
layer deepening, the correct final state is only achieved after ten days. The 1987 Mediter- 
ranean deep convection event was achieved over only a few hours. That was be the major 
difference found between the models and the in-situ data. With some adjustments in the 
tuning coefficients, the models can approach the rapidity of the actual event. But those 
models do not have the physics of the deep convection built in so the change of the pa- 
rameter's value may be irrelevant. 

We can recognize the ability for these simple models to produce believable long 
term results and need to look at 3-dimensional models. Among them, a non-hydrostatic 
large eddy simulation (LES) model (Moeng 1984) was adapted to ocean mixing. In the at- 
mosphere as well as in the ocean, the LES mode] attempts to simulate the large scale 
structure of the convection process (simulation of large eddies produced by the instabili- 
ties of the mean flow) and parameterize the smaller eddies that derive from bigger eddies 
by an energy-cascade process. Thus, solving the primitive equations (momentum, continu- 
ity, temperature and salinity) for the dynamic variables separated into large eddy fields 
(mean quantities) and subgrid-scale fluctuations (turbulent quantities), the LES model pro- 
vides the evolution of the velocity, temperature, and salinity fields as well as the statistic of 
the turbulent quantities. The first runs of this model under a constant forcing seem 


promising. 
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Figure 11. Kraus and Tumer Simulation of the Variation of Mixed Layer Depth, Tempera- 
ture, Salinity and Density with Time under Constant Windstress and Heat Flux. 
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Figure 12. Kraus and Turner Simulation of the Variation of Mixed Layer Depth, Tempera- 
ture, Salinity and Density with Time under Actual Windstress and Heat Flux. 
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Figure 13. Kraus and Turner Simulation of the Variations of Mixed Layer Depth, Tem- 
perature, Salinity and Density with Time under Actual Windstress (--) or Heat Flux (-.-). 
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Figure 14. Kraus and Turner Simulation of the Variation of Mixed Layer Depth, Tempera- 
ture, Salinity and Density with Time under Actual Windstress and Heat Flux using Ad- 
justed Coefficients (cl=1, c2=2 (....)/cl=2, Cm (-.-.-) /cl=2, c2#2 (----) ). 
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Figure 15. Naval Postgraduate School Simulation of the Variation of Mixed Layer Depth, 
Temperature, Salinity and Density with Time under a Constant Windstress and Heat Flux. 
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Figure 16, Naval Postgraduate School Simulation of the Variation of Mixed Layer Depth, 


Temperature, Salinity and Density with Time under Actual Windstress and Heat Flux. 
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Figure 17. Naval Postgraduate School Simulation of the Variation of Mixed Layer Depth, 
Temperature, Salinity and Density with Time under Actual Forcing and No Rotation. 
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Figure 18. Search for the Best Set of Coefficients. Naval Postgraduate School Simulation 
of the Variation of Mixed Layer Depth and Temperature with Time under Actual Wind- 
stress and Heat Flux using Different Sets of Coefficients. 
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Figure 19. Naval Postgraduate Schoo! Simulation of the Variation of Mixed Layer Depth, 
Temperature, Salinity and Density with Time under Actual Forcing using Different Coeffi- 
cents: ((M1, M2, M3, M4, MS) = (1, 0.5, 6, 1, 1)(----) / (0.5, 1, 12, 0.5, 0.5)(.....) / 
(0.25, 2, 24, 0.25, 0.25) (---) / (0.125, 4, 48, 0.125, 0.125) (-.-.-)). 
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Figure 20, Time Series of Mean In-situ Turbulent Kinetic Energy at 450 meters and Tur- 
bulent Kinetic Energy Produced by the Naval Postgraduate School Model under Actual 
Forcing using Different Coefficients (labelled as in Fig. 19). 


IV. REMOTE SENSING AND DEEP CONVECTION 


A. MOTIVATION 

EE Previous work on the subject 

During the past decade, satellites have proven to be very useful in measuring and 
understanding characteristic features of the ocean. Currents, sea surface temperature, sur- 
face winds, sea surface level, waves and biological activity are some examples where our 
knowledge has considerably improved through the use of satellite data. Altimeter, Syn- 
thetic Aperture Radar (SAR) and Advanced Very High Resolution Radiometer (AVHRR) 
are the most significant instruments. 

Detection of deep convection from orbiting platforms has rarely been attempted. 
European Remote Sensing Satellite ERS-1 SAR images of the Greenland sea have been 
compared (Carsey and Garwood 1993) to the Large Eddy Simulation (LES) model of 
Moeng (1984) adapted to ocean convection (Garwood et al. 1994) for the winter of 1992. 
The ice covered convection area (20 x 90 km) shows a "ragged net" pattern that was 
simulated by the model. The convection zone contains small, warm cells (~ 1 kilometer in 
diameter) of upwelling water surrounded by sinking, cold water. The cells rotate anticy- 
clonically and, in fact, are visible through the imprint they leave on the first-year ice. 
These encouraging results open the door to further investigation of deep convection proc- 
esses using satellites. 

To date, deep convection in the Golfe du Lion has not been examined using re- 
mote sensing. However, GEOSAT altimeter data are available for that area and studies of 
general Mediterranean sea surface heights have been conducted (Arabelos and Tziavos 
1990). 

Sea surface height studies have been conducted in all parts of the world ocean. 
Their results cover mainly currents (Gulf Stream), circulation (Southern ocean), bathymet- 
ric features (Gulf of Mexico) and oceanic events (El Nino). For more information, the 
reader should refer to the two Journal of Geophysical Research special issues (March and 
October 1990) devoted to GEOSAT altimeter results. 
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A attempt was made to look at AVHRR images of the Golfe du Lion during the 
convection period. Due to the cloud coverage only two images were available. The preci- 
sion of the radiometer as well as the software used to exploit the images were not suffi- 
cient to detect the temperature evolution of 1/10 degree related to the convection process. 

The investigation of deep convection using remote sensing instruments will be 
done with the GEOSAT altimeter data. 

2. Theoretical computation of the observable signal 

A preliminary and simple computation of the theoretical sea surface height during a 
convective event can be done using first a geostrophic computation and a mass conserva- 
tion type computation. The geostrophic (4.1) and hydrostatic (4.2) equations can be com- 
bined (4.3) to approximate the slope of the sea (n) surface level as a function of the 


surface velocity (u,v). 


=19P =-19P 
P= Bax Pay pee 
aP =pgan (4.2) 
ny AX 
foaeMn eA or an=f (4.3) 


With a surface circulation of 0.1 to 0.2 m/s and an elliptical shape area of deep 
convection whose major axis is 100 kilometers, the expected sea surface height variation 
should be 10 to 20 centimeters. 


The sea surface variation inferred from the mass computation is given by: 
XAtAG + 1)- AW] * [PGF 1) + P|DV2} = Paha (4.4) 


where h(i) [p()] is the depth [density] of the i® layer, h,, [p,,] the depth [density] reached 


after the convection (here bottom = 2000 meters). Using the density profiles before and 
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after the convection process (Leaman and Schott 1991) the sea surface variation should be 
(15 centimeters). These two simple computation show that the expected signal is at least 


greater than the altimeter accuracy, and therefore should be observable. 


B. THE GEOSAT ALTIMETER 

1. The geodesy satellite 

The GEOSAT (Geodesy satellite) was launched on March 12, 1985, into a 
108-degree inclination, 800-kilometer altitude orbit. The primary mission provided 18 
months of altimeter data required to improve the accuracy for the measurement of the 
earth's gravitational field. In October 1986, GEOSAT was maneuvered into a 17-day ex- 
act repeat orbit in order to detect mesoscale oceanographic features. The system then pro- 
vided a ground-track separation of 164 kilometers at the equator and 110 kilometers at 
midlatitudes (John Hopkins APL 1987). The GEOSAT program ended in October 1989 
due to a failure in the recording system. 

2. The principle of the altimeter measurement 

The principle of satellite altimetry is very simple. Short radar pulses are reflected 
from the Earth, The satellite height is deduced from the two-way travel time of the pulse. 
The altimeter measures the height of the satellite above the sea surface (Hs). 

‘The topography (Ht) is the departure of the ocean surface from the geoid and is in- 
fluenced by oceanic features, The orbit height (Ho) and the geoid height (Hg) are the 
heights of the satellite orbit and the geoid above an arbitrary reference surface, generally 
chosen to be an ellipsoid. To measure the ocean topography, we must combine the 
heights, 

Ht = Ho - Hs - Hg (4.4) 


The geoid and the time independent errors are a permanent signal (Hg) that can be 


computed by taking the mean of the satellite height (Ht+Hs) over a certain period. The 
orbit height (Ho) is obtained by a worldwide ground-based tracking network. 
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Figure 21. The Altimeter Measurement System. From John Hopkins APL, 1987 


3. The altimeter characteristics 

The GEOSAT altimeter is a 13.5-gigahertz nadir looking pulse compression radar. 
Estimation of three parameters is possible from processing the pulse compression data: 
height (distance between the satellite and the sea) with a precision of 3.5 centimeters for a 
2 meters significant wave height (SWH); SWH, with generally a precision of 0.5 meters 


and wind speed with 1.8 meter per second accuracy over the range | to 18 meters per 
second. 


4. Typical altimeter problems 

Instrumental and environmental corrections due to atmospheric water vapor, dry 
gases, ionosphere and ocean waves and tides need to be known to have accurate height 
measurements. Moreover, since the orbit is the reference frame (Stewart 1985) for the 


measurements, a special attention must be paid to the orbit error correction in altimeter 


computations. 





The orbit error is the sum of the radial error (error on satellite height), the 
along-track error (satellite ahead or behind its theoretical position on the orbit; this can 
also be treated as a height error) and the ground-track error (satellite left or right of its 
theoretical orbit). The clock error (sampling at an inappropriate time) can also occur but 
will be treated as an along-track error (Stewart 1985). Consequently the satellite orbit er- 
ror reduces to a ground-track and a height error. 

The ground track error can be solved using a crossover adjustment or an interpola- 


tion whereas the orbital error can be minimized by either crossover adjustment or a collin- 


ear technique. In all the cases, these types of error are difficult to model and to correct. 


The overall altimeter precision is the sum of all its errors and can be estimated to ten 


centimeters. 


Cc. DATA ANALYSIS 

1. The altimeter data 

The data were extracted from CDs-ROM available at the Oceanography Depart- 
ment of the Naval Postgraduate School. The CDs cover the entire Exact Repeat Mission 
(ERM): November 1986 through October 1989. Using the GEOSAT daily file, equator 
crossing list, orbit corrections file and mean tracks files as input, a FORTRAN program 
returns the height, time, latitude and longitude of the observations along each elementary 
track crossing the chosen domain. An option in the program enables the observed heights 
to be corrected for the wet and dry tropospheric error, the ionospheric error, the solid and 
ocean tides and the inverse barometer effect. Orbit correction does not exist and special 
attention must be given to this problem. 

2. The corrections applied to the data 

The orbit correction was estimated in two steps. First, a two year mean (Novem- 
ber 1986 - October 1988) of the global orbit was removed. Then the residual bias on each 
individual track was removed by subtracting the residual mean at that particular time and 


location. The importance of bias removal can be seen on Fig. 24 and 25. After this orbit 





correction, the remaining signal mainly lies between -10 and 10 centimeters with a few ex- 
ceptions consistent with the theoretical computation. 

In areas where the geoid gradient is large (such as the conelike topography of the 
Golfe du Lion), we must look carefully at the cross-track error. An easy way to check the 
size of this error is to compute the geoid with the 2 year mean tracks used for orbit cor- 
rection. This was done using the multiquadric interpolation program described later and 
produce a Western Mediterranean geoid close to the observed one. Therefore, we can 
conclude that the cross-track errors are small during this computation. 

Two atmospheric corrections were applied. The wet tropospheric correction was 
deduced from the measurements by the NOAA TIROS (Television and Infrared Observa- 
tion Satellite) operational vertical sounder (TOVS) or the Special Sensor Microwave Im- 
ager (SSM/I) on board the Defense Meteorological Satellite (DMSP). This wet 
tropospheric correction is induced by the water vapor present along the altimeter path. 
The propagation velocity of a radio pulse is also slowed by the dry gases present in the 
Earth's atmosphere. The chosen dry tropospheric correction is determined from a model 
assimilating data from the European Center for Medium Range Weather Forecasting 
(ECMWF). 

The inverse barometer effect, which is the hydrostatic response of the sea surface 
to the atmospheric pressure was assumed to be a constant over the small area considered. 

No tidal correction was made because no good tidal model was available for the 
Mediterranean. However, the error can be estimated. The M2 tide seems to be dominant 
in western Mediterranean sea and has an amplitude of 5-6 centimeters and a constant 
phase of 230° across the basin (Sanchez et al. 1992). 

Finally, even with all theses corrections, we must remember to be very careful 
when interpreting satellite data because of all unresolved or insufficiently modeled errors. 

3. Delimitation of the study and plan of investigation 

A long file of residual heights along with its respective time, location and track 


number is now available for further studies. Figure 22 shows the satellite ground tracks in 


the investigated area, With such data, the investigation of short time scale processes was 
impossible. The exact repeat of the satellite provided the same ground track only every 
seventeen days and the data were missing during certain periods. Due to the sampling rate, 
the distance between each observation is greater than 5 kilometers, so the altimeter sam- 


pling scale is greater than the small scale plumes of the deep convection. 
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Figure 22. The Area of Investigation and the GEOSAT Ground Tracks Used. 


However, a monthly study of the global variation of the sea level in a deep convec- 
tion area is possible. The preconditioning and sinking and spreading phase last a few 
weeks and can be associated with a sea level change. We will try to look at the large scale 
(O(50-100) kilometers) and long term (1 month) process occurring along with or because 


of the deep convection. 
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The computation was performed using the following strategy: 

© Regroup the data on a monthly basis. 

© For each month, compute mean tracks. 

© Run an interpolation program in between each mean track, and produce a 
monthly interpolated height field. 

© Produce a movie using the 24 interpolated frames. 

4. Producing monthly interpolated frames 

a. Monthly mean frames 

The first step of the computation was to regroup the data on a monthly ba- 
sis. The time of the sampling, given in seconds and referenced to January 1, 1985 was 
transformed into the Julian day. In order to simplify the computation, a month was 
rounded to 30.4 days. Then each line of the data was identified and placed into its respec- 
tive monthly frame. This produced twenty four frames (November 1986 to October 1988), 
each frame containing 5 to 16 tracks. 

Within each frame, the collinear tracks have been regrouped by identifying 
the track number available for each sampling. For each track number, all the data were 
read. Often the track was recorded twice, so it was necessary to remove duplicated data. 
If two similar times were found, one of them was discarded. The remaining collinear data 
were averaged. 

To compute a mean height along each collinear track, latitude is used as 
the reference axis, each track was divided into 1/10 degree bins (Fig. 23).Within each bin 
the latitude, longitude and height were averaged to produce one mean track (mean lati- 
tude, longitude and height every 1/10 degree). 

For each track, from one frame to another, the mean latitude and the mean 
longitude vary generally less than 0.05 degree (approximately 4 kilometers). Thus, it is 
reasonable to hypothesize that over the frames, each track has an identical ground signa- 
ture and therefore confirm that the cross track error is minimal in regions where the geoid 
has little gradient. 
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Figure 23. The Principle of Averaging Two Collinear Tracks. 
This computation produced 24 frames each containing 5 to 10 tracks with 
averaged residual heights. 


b. The interpolation process 

The field is sampled only at particular locations. We need an interpolation 
program to cover the entire domain. A multiquadric interpolation program developed 
Chumbinho (1993) following a method provided by Nuss (1993) was available. 

For each sampled location, the ground track (latitude, longitude) and the 
value of the surface (H, residual mean height) is known. The interpolation of the surface in 


any arbitrary point X is based on a weighted sum of radial basis functions: 
N 
H(X)= DaiQ(X-X,) (4.5) 
a 


, being the weight of the kernel function Q, and its argument the radial distance between 


the point X and the observation point X,. 
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For multiquadric interpolation, the kernel function is a hyperboloid 


function: 
Q(X -X)=[X-Xi)? +07}? (4.6) 


The constant c, an arbitrary small number, precludes the function from vanishing at the ob- 
servation point. 
Equation 4.5 is applied to all N observation points to solve for the weights 


producing a set of N equation with N unknowns. 
N 
Hy=H(X)) = LaiQilx,y)) (4.7) 
A 


‘The Q matrix is inverted to yield the weights that in turn provide the inter- 
polated solution at any regularly spaced grid point: 


Hy = QgiGt; where Ogi = [(xg— xi)? + gyi)? +07}? (4.8) 


where (x, ,y,) are the coordinates of each regular grid point. The solution depends only on 
the number of observation points available and is independent of the size of the grid to be 
interpolated. 

To account for noisy observations an uncertainty parameter and a smooth- 
ing coefficient are introduced in the diagonal element of the kernel function (Nuss and 


Title 1993). 
Hy =Qy+[N*A* 07 * 84) * a; (4.9) 


where 8, is the Kroneker delta function, 0,” is the smoothing coefficient and 2 is a vector 


of uncertainties in the observations (mean-squared observation error). 
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c. Searching the best interpolation coefficients 

Choice of the best smoothing parameter and coefficient c is somewhat sub- 
jective. The damping of the field is controlled by the smoothing parameter; large values 
cause a strongly damped solution. The constant parameter determines the shape of the in- 
terpolated field; large values of c imply a rather flat interpolating function, and tight gradi- 
ents or closely spaced observation points may be misrepresented. The table below 


summarizes the different choices applied in conjunction with Fig. 26. 

















{LARGE LAMBDA SMALL LAMBDA 
Field strongly damped, —_[ Real field. 

LARGE C Misrepresentation of Misrepresentation of 
a tight gradient. a tight gradient. 
Field strongly damped. _|Real field. 

SMALL C Good representation of _ |Good representation of 
a tight gradient. a tight gradient. 














Table 2. Influence of the Interpolation Coefficients on the Field Produced 


Inappropriate values of the smoothing parameter create an unrealistic inter- 
polated field. The initial heights show many bumps along their tracks that need to be 
smoothed. However, the general shape of the field must be preserved without too much 
damping or contamination by false data created during the interpolation (black spots on 
the plots are points out of range). As it can be seen on Fig. 26, using different pairs of pa- 
rameters leads to different interpolation fields. The pair chosen for the study is the best 
trade-off possible with a medium constant parameter (c=10*) and a small smoothing coef- 


ficient (A=10"). This allows reasonable representation of the general gradient. 
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D. RESULTS 

For each of the two runs, the full area and a zoom over the area of interest (100 
kilometers around the center of deep convection: 42°N - 5°E) were presented. In the first, 
the averaged heights were used without smoothing. In the second case, for every month, a 
mean shape was visually extracted from each individual track. This was done to reduce the 
bumps found on each tracks, because the sea surface rarely has steep gradients. The plots 
are presented at the end of this chapter. 

1. The averaged heights 

a. Interpretation of the results 

Looking at the global area (north part of the Northwestern basin), there is 
unfortunately no direct correlation apparent between the deep convection and the evolu- 
tion of the sea surface level (Fig. 27, 28, 29 and 30). Except for a small feature in the 
northwest of the first three frames, there is no obvious common feature that could be seen 
when playing the movie made with the monthly frames. There is also no direct comparison 
from one year to the next. The rise and fall of the sea level neither occur at the same loca- 
tion nor at the same time in the year. 

However, this small feature deserves a special attention (Fig. 30). During 
the three months prior to the convection period (November 1986, December 1986 and 
January 1987) there is a significant deepening close to the estimated location of the con- 
vection (42°N - 5°E) which happens to be along satellite track 487. The elevation of the 
convection area (100 kilometers around that location) is lower from November to January 
and appears slightly higher afterwards. This late winter deepening may be linked with the 
preconditioning phase and the early spring rising with the spreading phase, when the water 
rises slowly from the bottom. 

Looking at that particular track (Fig. 31), the presence of the low can be 
seen in almost all of the passes and seems to be moving between 41.5°N and 42°N. Since it 
appears on the original tracks, the feature is not a series of false data created by the pro- 


gram, so it can be interpreted as a real phenomenon. 
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b. Comparison to in-situ observations 
How does the above discussion relate to the investigation of Schott and 
Leaman (1991)? The position of their mooring was based on the supposed location of the 
center of the deep convection; they also conducted CTD observations up to 60 kilometers 
from the center to determine the hydrographic structure within the convection zone. Their 
data cover only the convective period (mid January 1987 to mid March 1987) and the 
closest satellite tracks only sample the southern portion of their zone of investigation 
(tracks 487 and 160). A direct comparison for the global area is difficult; however, the po- 
sition of the feature discussed earlier coincides with the southern limit of the deep-mixed 
patch (salinity > 38.4) found by Schott and Leaman (Fig. 32). Part of the signature of this 
heavier water could be the low seen in the satellite interpolated field. 
2. The smoothed tracks 
To reduce the bumps in the sea surface heights, a mean shape was visually ex- 
tracted from each individual track. This was done to try to reveal the large scale feature. 
The interpolated field built with these smoothed tracks is flatter than before (Fig. 33) and 
highs and lows are now easily identifiable. As with the previous field, it is not possible to 
firmly associate the evolution of the convective patch with some of the images. The previ- 
ous interesting feature shows up only in January 1987 but it is now larger. The smoothed 
tracks do not bring any additional information. 
3. The averaged winter frame 
A final and probably better use of these frames can be made by averaging all the 
winter frames and comparing them to the mean winter surface density field. From the 
analysis of historical data (1905 to 1988) a winter sea surface density field was deduced 
(F. Aikman III, personal communication) which helps to define the extend, as seen in 
Fig. 34, of the deep water formation in the Northwestem Mediterranean sea (Madec et al. 
1991). 
Frames from November to March were used to produce a mean sea level field 


that shows a characteristic slope toward the north (Fig. 35). The blue oval shape area 
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refers to a mean low that is fairly consistent with the limits of the heavier water. As we go 
south, the color changes from blue to yellow, meaning that the sea level increased. This 
also means that the water is becoming less dense. The blue area is slightly higher in lati- 
tude than the historical data field but our computation uses only two years of data instead 


of 83 and there are still unresolved errors. 


E PARTIAL CONCLUSION AND RECOMMENDATIONS 

For the large area and the two year period covered, it is not possible to extract the 
location and the extent of the deep convection, nor firmly identify it in time. However, 
knowing that there was a deep convection event, there is a small feature found over three 
months in the northwest of the remotely sensed area which could be related to the low in 
the sea surface. A correlation with in-situ data is possible. 

The sum of the residual errors and the insufficient spatial and temporal resolution 
prevent the observation of the typical deep convection lows. The tidal and inverse barome- 
ter errors which are not taken into account and the orbital error which is difficult to esti- 
mate are certainly major handicaps. 

We can go further in the investigation only if we possess the correct tools in terms 
of satellite instruments. An increase the spatial and temporal resolution are strongly re- 
quired in order to see small-scale and short-period phenomena that might be identify as 
deep convection. The TOPEX/POSEIDON altimeters will probably give a significant im- 


provement to the remote sensing of oceanic convection and related mesoscale features. 
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Residual heights after a 2 year mean global orbit correction 
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Figure 24. The Residual Heights after the Two Year Mean and Bias Removal. 
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Track 487 after a 2 year mean global orbit correction 
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Figure 25. The Residual Height of Track 487 after the Two Year Mean and Bias Removal. 
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Figure 26. The Influence of The Multiquadric Interpolation Parameters on the Interpolated 


Field. 
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Figure 27. Monthly Frames (November 1986 - July 1987) of Interpolated Field Produced 


with Averaged Heights (Scale in Meter). 
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Figure 28. Monthly Frames (August 1987 - April 1988) of Interpolated Field Produced 
with Averaged Heights. 
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Figure 29. Zoom of the Convection Area (November 1986 - July 1987) of Interpolated 
Field Produced with Averaged Heights. 
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Figure 30. Same as Fig, 29 but for the Period November 1986 - February 1987. 
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Figure 31. Time Series of the Track 487 (November 1986 - February 1987). 
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Figure 32. The Deep Mixed Path inside the Convection Area during the Convection Pe- 


riod. From Leaman and Schott 1991. 
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Figure 33. Monthly Frames (November 1986 - July 1987) of Interpolated Field Produced 
with Smoothed Heights. 
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Figure 34. Sea Surface Densit 
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Figure 35. Mean Winter Interpolated Altimeter Field. 
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Vv. CONCLUSION 


The first part of this research investigated the behavior of the Kraus and Turner 
and the Naval Postgraduate school mixed layer models under deep convection conditions. 
These models were initialized with observed temperature and salinity profiles and were run 
under the observed forcing. With the tuning coefficients set for a nonconvective mixed 
layer deepening (original set of coefficients), both models produce final values of tempera- 
ture and salinity similar to observations. However, the mixed layer depth only reached 
1150 meters with the Kraus and Turner model and 500 meters with the Naval Postgradu- 
ate School model after ten days of forcing. The swiftness of the deep convection process 
(a few hours in the observed case) can be approached by adjusting the coefficients. In that 
case, the Naval Postgraduate School model provides the best simulation and reaches the 
bottom and the observed temperature and salinity deep convection values after three days. 

The changes in the coefficients are required by the convection process, but a physi- 
cal justification can be made for only two of them. A complete justification of the coeffi- 
cients can be done with higher order convection models and the Large Eddy Simulation 
model may be used to better understand the physics of deep convection. 

An investigation of deep convection by remote sensing instruments has been con- 
ducted in the second part of this research. The identification of this process with an inter- 
polated altimeter field is still not possible on a global temporal and spatial scale. There is 
no strong signature of the deep convection process. The required signal magnitude is 
slightly higher than the precision of the measurement, and some errors are still present. 
However, a persistent low in the sea level at the location of the deep convection area and 
in the mean winter altimeter field are consistent with in-situ observations and provide 
some hope in studying and later identifying this phenomenon from space. With the recent 
altimeters, such as those on the joint United-States/France TOPEX/POSEIDON mission, 


we may now have a tool to answer this investigation. 
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